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Summary 


NASA Glenn Research Center has developed MathWorks® 
Simulink® controller logic for driving a room-temperature, 36- 
tooth stator, four-pole rotor reluctance motor. The Simulink® 
core logic was extracted from an existing C motor controller 
that was employed to achieve a rotor speed of 3000 rpm. The 
Simulink® controller, however, has significant additional logic 
refinements that were not available in the existing C controller, 
such as the “Per Rev Logic” component and its frequency 
filter. The filter provides a more accurate reading of the rotor 
input signals than the previous version of the controller. The 
controller is versatile, and with slight modifications, can be 
used to drive other reluctance motor types incorporating 
dissimilar stator and rotor pole combinations. The original C 
controller was designed with the goal (after appropriate 
modification) of controlling a future superconducting motor. 
This superconducting motor will be employed as a testbed for 
developing other superconducting aviation propulsion motors 
envisioned for future turboelectric aircraft. The Simulink® 
results presented in this paper were generated from simulated 
rotor inputs. However, in an actual application, these simulated 
inputs would be replaced by actual proximity probe signals. 


Introduction 


A controller was designed for a motor that will eventually 
Operate as a wound rotor synchronous machine. The rotor has 
four salient poles on which superconducting coils will 
eventually be installed (Fig. 1). However, it was initially 
operated without rotor coils (Fig. 2) as a reluctance machine. 

The 36 stator teeth (Fig. 3) are each wound with a 
concentrated coil that encircles only one tooth. Either 
configuration of the rotor displayed above can be inserted into 
the stator. 

For reluctance mode operation, the coils are connected in 
nine separate phases, each with an independent power 
converter. However, for the simulation discussed in this paper, 
only three phases were utilized, so only every third coil 1s 
activated. Each phase has four coils in series, equally spaced 
around the stator, with one coil in each quadrant of the stator. 
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The width of a rotor pole is sufficient to overlap no more than 
three stator teeth at a time, only one of which 1s active in this 
study. So at any instant, many of the potentially active stator 
teeth are not near any rotor pole. Energizing coils, when they 
are too far from a rotor pole, would waste power without 
producing torque. Consequently, one aim of the controller is to 
send current to a coil only when a rotor pole is approaching and 
to turn the current off before the trailing edge of the rotor starts 
to pass the stator tooth and generate torque in a backward 
direction. Strong torque is generated by any stator tooth only 
when it is energized at a time when a rotor pole is approaching 
or partially overlapping the tooth. Therefore, the highest 
average torque per dissipated power is expected to be produced 
with current waveforms that are more pulse-like than sinusoidal. 

Thus, from the configuration of the motor and the nature of 
its torque production, the controller must be able to command a 
pulse-like current waveform for each stator phase at a time 
when the four rotor poles are approaching and overlapping the 
four coils in the phase. Hence, determining rotor position is the 
first task of the controller. Based on the instantaneous rotor 
position, the controller should send the appropriate commands 
for energizing each of the phases. Most of the phases will not 
be significantly excited (1.e., will be “off’) at any instant. The 
electrical control commands from the dSPACE GmbH 
hardware are sent to three sets of four stator coils. 

Figure 4 shows a block diagram of the system. The four coils 
in a set are connected in series to produce magnetic pole 
combinations north-south-north-south (N-S-N-S) and are 
spaced 90° apart in the stator. In a three-phase configuration, 
each coil set 1s separated 120 mechanical degrees from the 
other. The coil combinations in the three phases are {3, 12, 21, 
303, {6, 15, 24, 33}, and {9, 18, 27, 36}. These coil sets together 
produce torque on the four-pole rotor. 

The rotating rotor produces electrical pulses from two 
proximity probes secured on the motor containment structure. 
The sensor end of one probe is positioned to sense the passage 
of two equally spaced notches on the shaft. The second probe 
senses the passage of each tooth of a 16-tooth gear mounted on 
the shaft. The combination of the pulses from these probes 1s 
sufficient for accurate timing of control commands for the stator 
coil power amplifiers. 


Figure 1.—Four-pole rotor with superconducting coils. 


Figure 2.—Four-pole rotor without superconducting coils. 
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Figure 3.—36 stator teeth with concentrated coils. 
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Abbreviations and Acronyms 


AATT Advanced Air Transportation Technologies 

GUI graphical user interface 

PAC Phi Angle Controller 

PID Proportional Integral Derivative 

PRL Per Rev Logic 

RPM ~ Revolutions Per Minute or revolutions per minute 
WSG Wave Shape Generator 


Symbols 
a variable, where a= —_ 
2 
B magnetic field flux density 
C pulse width 
do differential angular displacement 
dV differential flux volume 
dw differential magnetic energy 
g gap length 
I rotor tooth width 
i stator coil current 
L period 
is lamination stack length 
n number of iterations 
N number of turns in the stator coil 
Dp sine function exponent value 
r outer radius of rotor tooth 
T instantaneous torque 
t time in seconds 
V flux volume 
Ww magnetic energy 
x the argument in a trigonometric function 
0 angle 
Lo permeability of free space 


PWM 


amplifier 


Current output 


°|Command 


Motor 


input 0 lL. One 
stator | ; 
four-coil }~ = 1 ofrine 
Control & wee —_— ! | Gees 
computer bt Current cS}! | group 
sense ~ | | 
board —)|! | 
CC a | 
! 


Not shown are liquid nitrogen fill controls, 
thermocouples, and one-per-rev sense circuit. 


Figure 4.—High-temperature superconducting motor wiring. 
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Rotor-Stator Mechanical Depiction 


Figures 5 to 8 show a combination of room temperature and 
superconducting rotor assemblies. The controller discussed in 
this paper was initially designed for the room temperature 
version of the motor, with the intent of expanding the logic to 
operate a superconducting motor configuration. Figures 7 and 8 
depict the rotor configuration implemented for the original C 
controller. The original C code that has been converted to its 
MathWorks® Simulink® version will be discussed in detail. 

Figure 9 depicts the final prototype assembly. The signals in 
the stator coils produce torque on the rotor’s four poles 
(see Fig. 2), thus rotating the rotor. The motor’s behavior 
approximates that of a synchronous motor, with the rotor 
rotation dutifully following the frequency of the sinusoidal 
signal in the stator poles. The instantaneous speed of the rotor 
is continuously compared to the input desired speed in the “PID 
Controller” Simulink® block, thus effectively maintaining the 
rotor torque and speed at a desired input value, and hence 
closing the control loop. The theoretical torque generated at a 
rotor-stator pole pair is derived in the General Magnetic Torque 
Equation section. It is an approximation that ignores the torque 
contribution of the fringing field, but is sufficient to allow 
preliminary evaluation of the effectiveness of a candidate 
waveform for stator coil excitation. 
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Figure 5.—Exploded view of 
superconductor rotor assembly. 
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Figure 7—Rotor fully assembled. _ 


Figure 9.—Prototype assembly configuration. (a) After rotor 
insertion (external view). (b) Bottom bearings with flange 
installed after rotor insertion. 


General Magnetic Torque Equation 


Figure 10 depicts an approximation of a portion of the rotor 
and a portion of a stator tooth on which the stator coil is 
energized to produce a magnetic field flux density B. The 
magnetic field interacts with the rotor at a gap length g. 
Differential flux volume is represented by dV. The magnetic 
energy, w, 1S given by the volume integral equation (Ref. 1). 


w= 


= | | B2dV (1) 


where L, 1s the permeability of free space. Integration over only 
the flux volume V is possible if the magnetic core material is 
operating within its linear range so that the energy stored in the 
core can be neglected compared with the energy stored in the 
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gap. Referring to Figure 10, it can be shown from Equation (1) 
that the instantaneous torque on the rotor, 7, towards the stator 
is given by 


2 
— Lsgr increasing overlap volume 
2 Lo 
dw 
Tf =-—= 0 constant overlap volume 
do ; (2) 
Lsgr decreasing overlap volume 
20 
where 


dw is the differential magnetic energy 

d@ is the differential angular displacement 

L, 1s the lamination stack length going into the paper 
g is the gap length 

r is the outer radius of the rotor tooth 


Furthermore, 
T = (uoN*Lsgr ji (3) 


where WN is the number of turns in the stator coil and i is the 
current in the stator coil. 


Fringing B 
field lines—~_ 


Figure 10.—Side view sketch of rotor-stator segment. B, 
magnetic field flux density; g, gap length; /, rotor tooth width; 
r, outer radius of the rotor tooth; dV, differential flux volume; 
dé, differential angular displacement. 
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Figure 11.—Torque potential (blue) on rotor superimposed 
on adjacent stator teeth coil current (red). 


In previous tests, it was observed that the rotor responded 
favorably (torque and smoothness of rotor revolutions per 
minute (rpm)) to drive signals that are representative of sine 
signals of odd power exponents (see the red curve in Fig. 11, 
Sine (8)9, or triangular wave (Fig. 13(j)). The torque potential 
(blue), and coil current (red) curves versus nondimensional time 
are depicted in Figure 11. Referring to Figures 10 and 11, the 
maximum clockwise torque on a rotor’s pole 1s achieved when 
the apex of the current pulse (red) in the stator coil 1s timed to 
engage when the leading edge of the rotor pole is positioned just 
at the trailing edge of a stator tooth, and vice versa. In addition, 
the width of the root of the current pulse should be such that its 
leading edge returns to zero value as the rotor trailing edge 
comes into alignment with the stator trailing edge. If the coil 
has any residual current at this instant, it will have an unwanted 
braking effect on the rotor due to the presence of fringing B 
fields that produce torque in the reverse direction. 

The coil current remains at zero as the subsequent rotor pole 
approaches the stator teeth, whereby the coil current then 
reverses its direction and produces torque on the rotor as 
described above. The current reversal switches the magnetic 
polarity in the stator teeth, thus having the beneficial effect of 
demagnetizing stator teeth and simultaneously maintaining the 
induced magnetic polarity in the rotor tooth. This electrical 
event occurs in all four coils in a set, and all three coil sets are 
energized simultaneously, thus producing continuous torque on 
the rotor. Table I depicts the calculated benefits (with the help 
of Mathematica) of using sine drive signals of odd exponents. 
The merit value is obtained by dividing the average torque by 
the expended energy. 

The expended energy does not consider eddy-current losses 
in the iron core, which increase with frequency and would make 


very sharp pulses less desirable. Incorporating those losses is 
outside the scope of this paper. Notice that at Sine (8)p, p = 9, 
the average torque is at its maximum value for a relatively small 
expenditure of energy, after which the torque begins to 
diminish. 


Graphical User Interface 


The graphical user interface (GUI) is depicted in Figure 12. 
For the rest of this report, texts that are enclosed in quotation 
marks represent the names of Simulink® blocks, inputs, and 
outputs as they appear in the GUI. The Simulink® controller 
presented in this paper uses simulated proximity probe inputs 
that can ultimately be replaced by actual proximity probe 
signals from dSPACE hardware. The actual rotating rotor 
produces regularly pulsing proximity probe signals that vary in 
frequency in proportion to the rotor’s speed. The pulses are 
created from 2 notches or 16 notches, located on the rotor’s 
shaft. The 2 notches are located 180° apart from each other, 
while the 16 notches are separated 22.5° apart from each other. 
Referring to Figure 12, the simulated input pulses are converted 
to “RPM,” and angle “THETA” (radians) in_ the 
“Per Rev Logic” block. The “RPM” signal is sent to the 
“Phi Angle Controller” block, and the “PID Controller’ block, 
while the “THETA” (radians) signal is forwarded to the 
“Wave Shape Generator” block for use as the argument in 
trigonometric functions. 

Efficient motor operation and torque are affected by the type 
of drive or control signal that 1s applied to the stator coils: 
“Phasel,” “Phase2,” and “Phase3,” via (SPACE hardware (see 
Fig. 12, midcenter of the GUI). Hence, 10 trigonometric drive 
functions or wave shapes are provided for testing and analysis to 
determine which signal gives optimal performance during 
motoring operation. The “PID Controller’ regulates the 
amplitude of the 10 trigonometric drive functions or signals 
residing in the “Wave Shape Generator” block. These functions 
are, namely; Sine, Sinu Pulse, Sine“, Trapezoidal 01, 
Trapezoidal_ 02, Square _Wave_01, Square _Wave_02, Sawtooth, 
Triangle 01, Triangle 02. These functions are italicized for the 
remainder of this report. 

Any of these functions can be selected to provide its 
characteristic waveform electrical signal, which is then 
outputted simultaneously from the “Phasel,” “Phase2,” and 
‘“Phase3” ports located on the “Wave Shape Generator” block. 
The phases are separated 120 electrical degrees. These output 
phase leads are coupled to the stator coils via dSPACE 
hardware. 


7 I. — VERSUS ENERGY EXPENDITURE FOR SINE DRIVE SIGNALS OF ODD EXPONENTS 


0.082 0.137 0.149 0.151 0.150 0.147 0.142 0.124 0.070 0.024 


0.0076 


Petey 0.314 | 0.169 | 0.155 | 0.132 | 0.117 | 0.106 | 0.091 | 0.066 | 0.031 0.011 0.0034 
expended 


“Sine function exponent value. 
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Pulser 


The simulations discussed 1n this paper employed input signals 
emanating from the signal generators “Pulserl,” “Pulser2,” and 
“0...6000 rpm” block (see Fig. 12, right center of GUI). The 
selection of any of these generators is facilitated by “Sw5,” 
“Sw6,” “Sw7,”’ and the “Pulser Selector.” The 
“Phi Angle Controller” and “PID Controller” were disabled 
(using “Switch” and “Sw4,” respectively), thus allowing 
observation (“Scopel”’) and optimization of the code settings, by 
way of the pure, zero-phased, and unmodulated signals (see 
Fig. 13) generated by the “Wave Shape Generator.” More 
information on the “Wave Shape Generator’ 1s available in the 
Wave Shape Generator section. 
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Referring to Figure 12, “Pulser 1” simulates 2 pulses per 
revolution (2 pulses per second is converted to 1 Hz of rotor 
rotation), while “Pulser 2” provides 16 pulses per rotor revolution 
(16 pulses per second are converted to 1 Hz of rotor rotation). 
Figure 14 depicts the subsystems of “Pulser1l” and “Pulser2” and 
depict the operational manifestation of these pulser signals. The 
pulsers generate only fixed frequencies, the values of which are 
entered in the MathWorks® MATLAB? interface as “f = #” Hz. 
In actuality, the rotating rotor produces regularly pulsing signals 
that vary in frequency with the rotor’s speed. This signal modality 
is more realistically achieved with the “0...6000 rpm” simulator 
signal block. The subsequent section discusses this quasi-realistic 
signal generator block when used in conjunction with the “PID 
Controller.” 
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Time, s 
Figure 13.—Wave_Shape_Generator signals. (a) Sine (8). (b) Sinu_Pulse. (c) Sine (8)‘p, p = 9. (d) Trapezoidal_01 
(e) Trapezoidal_0O2. (f) Square_Wave_07. (g) Square_Wave_02. (h) Sawtooth (2Hz). (i) Triangle_07. 


(j) Triangle_02. 
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Scope 


“4% 


Pulser 1 Sine Wave' > 
Math r 
Function3 a4 
For Sine Wave1 g 
Amplitude: 1 
Bias: 0 Constant 


Frequency (rad/s): 2 x pi x f 
Phase (rad): 0 
Sample time: 0 


(a) 


Sine Wave1 


Math 
Function2 


Voltage, v 


For Sine Wave1 
Amplitude: 1 
Bias: 0 
Frequency (rad/s): 2 x pi x 8 x f 
Phase (rad): 0 
Sample time: 0 


(b) 


Constant 


.O 
Re 


oO 
on) 


2 pulses per complete revolution 


0.0 0.1 


0.0 0.1 


02 03 04 0.5 06 0.7 08 0.9 1.0 
Time, s 


02 03 04 0.5 06 0.7 08 0.9 1.0 
Time, s 


Figure 14.—Pulser subsystems. (a) 2_Per_Rev_Pulser (1 cycle x 2). (bo) 16 Per _Rev_Pulser (8 cycles ~ 2). 


0...6000 rpm 


To simulate the rotor’s input, the “O...6000 rpm” block was 
used to mimic variable rotor speeds over a | s period. This 
artificial rotor speed vs. time is chosen to illustrate the resulting 
command output of the controller; 1t is not the speed history that 
would result from controller output. The simulation was done 
with the “PID Controller’ enabled via “Sw4” and the 
“Phi Angle Controller’ disabled (“Switch”). Note that during 
actual motor operation, the “Phi Angle Controller” is enabled 
(the enabled “Phi Angle Controller” at —100° phase shift 
voltage graph differs only slightly from the disabled 
“Phi Angle Controller” voltage diagram, and was excluded). 
Figure 15(a) illustrates an inputted 1 s duration simulated 
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variable rotor speed signal originating from the “0...6000 rpm” 
block. The block’s signal reached the 3000 rpm maximum limit 
and returned to 0 rpm. The corresponding control voltage signal 
is depicted in Figure 15(b). In Figure 15(b), the “PID 
Controller’ is observed continuously adjusting the voltage level 
of function 2—Sinu Pulse (see the “Fn Select” setting), as it 
attempts to limit the rotor speed to 500 rpm, which was the 
initial desired speed. Note that the simulation setting can also 
be changed to 6000 rpm by simply assigning “Amplitude = 
10.4414” in the MATLAB® interface or by using appropriate 
dSPACE settings (the “PID Controller” settings may need some 
adjustment). The subsystem of “O0...6000 rpm” is displayed in 
Figure 15(c). The “Per Rev Logic” block will be discussed 
next. 


0.0 0.1 0.2 0.3 04 05 06 0.7 08 0.9 1.0 0.0 0.1 0.2 0.3 04 0.5 06 0.7 08 O39 1.0 
(a) Time, $s (b) Time, s 


For 3000 rpm max, set amplitude in sine function to 7.05. 


2.983e-30 For 6000 rpm max, set amplitude in sine function to 10.4414. 


rpmset (rpm) For Sine Wave1 and 2 


Amplitude: Amplitude x 0.8765 
Bias: 0 
Frequency (rad/s): (2 x pi) 
Scope Phase (rad): 0 x (pi/180) 


revolutions / min Hz 2 + Sine Wave1 ramp 1 
; 
rpm 
Math Add 
revolutions/min Function 


Sine Wave 
For Sine Wave Function 1 and 2 Function’ 
Amplitude: 1 
Bias: 0 


Frequency (rad/s): (2 x pi x 4) 
Phase (rad): 0 
Sample time: 0 


Sine Wave2 


rads/sec ramp 2 


rads/sec 


Math Add1 


rads/sec Function 


Sine Wave 
(c) Function2 


Figure 15.—0...6000 rpm block. (a) 0 to 3000 rpm and back to 0 rpm. (b) Control voltage level for 0 to 3000 rpm 
and back to 0 rpm rotor speed. (c) 0...6000 rom subsystem. 


NASA/TP—2017-219451 11 


Per Rev Logic 


This section describes the features of the “Per Rev Logic” block. The “Per Rev Logic” subsystems’ level progression 1s 
shown in Figures 16 to 19. Figure 19 depicts the “RPM” blocks “2_per_rev” and “16_per_ rev.” These blocks contain the 
subsystem that converts the proximity probe pulses, generated by the rotor’s notches, to the corresponding rotor rotational 
frequency. The subsystems are structured in this order, from top to bottom: 


e “Per Rev Logic” 
o “Per Rev Logic 3” 
“Per Rev Logic 2” 
@¢ “Per Rev Logic 1” 
» 2 and 16 Per rev 


Per_Rev_Logic_3 


(Swit [Switch ED 


Pulser Selector 
From1 Goto1 


instantaneous (radians) 


THETA Pulser‘ 
: 
RPM Pulser2 
From2 
Multiport Per_rev_logic 3 
Switch1 
instantaneous (rpm) 
Figure 16.—Per_rev_logic subsystem. 
Per _Rev_Logic 2 
50.27 
radians1 pulsert 
; Selector 2 per rev 
7 
theta (rads) 
16_per_rev 
Multiport 
Switch 


Per_rev_logic_2 


radians2 


pulser2 


Figure 1/.—Per_rev_logic_3 subsystem. 
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Per Rev_Logic_1 
it 


Rads 1 
Ramp1 
rads1 
SLU dina Radians1 
rpm rad radians/sec_1 
Scope Scope2 
(2) 
RPM 1 sensor_signal 1 
oy 
RPM 2 Per_rev_logic_t sensor_signal 2 
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Figure 18.—Per_rev_logic_2 subsystem. 
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Figure 19.—Per_rev_logic_1 subsystem. 
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The “Frequency Filter1” block is located one level down in 
the “2 per rev” subsystem block (Fig. 20). This block 1s 
located at the top right corner of the frame. This filter removes 
intrinsic frequency anomalies from the drive signal. These 
anomalies occur prior to the “Divide” block when the “Digital 
Clock” and “Delay” block attempt to calculate the time 
differences between pulses at the initiation of signal detection. 
This anomalous behavior is also evident in the “16 per rev” 
subsystem (Fig. 21). These subsystems are triggered at each 
signal pulse. During the fractions of a second prior to the time 
difference (Af) calculation, the time value is significantly 
smaller than the At used to divide the “angle between_teeth” 
value. The division by this smaller time value results in an 
elevated (anomalous) value of degrees and seconds, and hence 
the frequency (see Fig. 22). This graph (Fig. 22) is clearly 
different from Figure 13(b). 


After the passage of several more fractions of a second, the 
correct At value is attained. Figure 22 should depict a 1-Hz 
signal; instead, there are 16 pulses in approximately 
0.5s. The accurate arithmetic at this juncture indicates a 
frequency of approximately 31.36 Hz. This anomaly persists 
across the entire frequency range, up to and including 100 Hz 
(6000 rpm), wherein the initial control pulse for each frequency 
value is anomalous. The anomalous values for a given set of 
predetermined frequency values ranging from 0 to 100 Hz were 
meticulously logged (see Table IT), and the values were graphed 
(see Fig. 23). It should be noted that for rotors with large mass 
and/or inertia such as the one discussed in this paper, this 
frequency anomaly does not affect motor operation. This is only 
discussed here to be thorough. 
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Figure 20.—2 Per Rev subsystem. Time is t. 
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Figure 21.—16 Per _Rev subsystem. Time is t. 
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Figure 22.—Wave anomaly for 1-Hz signal. 
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The Table II values were subsequently entered into a one- extrapolation method. Figures 24 and 25 are the subsystems of 
dimensional lookup table function (see “Filter 1,” Fig. 24; the “Frequency Filter 1” and “Frequency Filter 2” blocks, 
“Filter 2,” Fig. 25). The “Cubic spline” was selected as the respectively. 
interpolation method, and “Linear” was chosen as _ the 


TABLE II._FREQUENCY FILTER VALUES 
FOR 0 TO 3136 HZ (Hz_1) AND 
0 to 100 HZ (Hz._2) 
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Figure 23.—Frequency filter. 
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Figure 24.—Frequency Filter 1 subsystem. 
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Figure 25.—Frequency Filter 2 subsystem. 
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Phi Angle Controller 


With the filters activated, a simulated 2-pulse-per-second 
(60 rpm) rotor signal from “Pulserl” was introduced into the 
“Phi Angle Controller,” first with the “Phi Angle Controller” 
disabled (“Switch’ off) and then with _ the 
“Phi Angle Controller” enabled (“Switch” on). This action 
applied a simulated -100° phase shift (acceleration 
“Phi Angle’) adjustment to the control signal (see Fig. 26(a), 
without phase shift and Fig. 26(b), with phase shift). 

This angle could well have been determined to be +100°. The 
angle’s ultimate sign (— or +) and magnitude are determined 
during the initial slow speed experimental setup of the motor 
rig. The —100° angle was automatically applied by the code to 
accelerate the rotor from 0 rpm to the desired initial speed of 
500 rpm (see Fig. 12, “Input (rpm)” top center of GUI). Note 
that this angle could have been determined experimentally to be 
a deceleration angle for function 2 (see “Fn Select’) during the 


1.5 
1.0 
0.5 


0.0 


Voltage, v 


D 
0.0 0.1 


Voltage, v 


5 
0.0 0.1 


(b) 


initial slow speed motor operation setup. The “Phi Angle” 
value will often differ for the other drive signals. As the rotor 
speed approaches, then exceeds 500 rpm, the “Phi Angle” 
switches to 0° (a simulated deceleration angle). A similar 
process occurs for the other functions at a given rotor speed, 
wherein the “Phi Angle” value may be different. 

The “Phi Angle Controller’ enhances the functionality of 
the “PID Controller’ when both are simultaneously engaged. 
The “Phi Angle Controller’ continuously monitors the rotor 
speed “RPM” and modifies the electric phase angle of the 
control signals with a leading or lagging “Phi Angle” value 
(relative to the stator poles). When the rotor reaches the desired 
speed, the “Phi Angle Controller” adapts a semineutral stance 
until 1t senses the slightest change in rotor speed. Note that this 
action or reaction occurs in microseconds. During this 
semineutral stance, the “PID Controller” actively maintains the 
desired speed of the rotor until a new speed is entered, and the 
control action repeats. 


0.2 03 04 05 06 0.7 08 09 1.0 
Time, s 


02 03 04 05 06 07 08 09 1.0 
Time, s 


Figure 26.—Phi_Angle. (a) At 0° (without phase shift), for f= 1 Hz, 
Proportional_ Integral_ Derivative (PID) disengaged. (b) at —100° 
(with phase shift), for f= 1 Hz, PID disengaged. 
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For this motor, the Sine (8)p, p = 9, and Triangle 02 waves 
(see Figs. 13(c) and 13(j)) appear to give the smoothest 
operation throughout the speed range of 0 to 3000 rpm, which 
is the maximum upper design limit of the rotor. The controller 
is capped at 6000 rpm, the maximum speed of a future rotor. 


Figure 26 depicts the actions of the “Phi Angle Controller” 
with the “PID Controller” disengaged and using an amplitude 
of 1 V (see Fig. 12, “Sw4’’) and using function 2 (“Fn Select’’). 
The subsystems for the “Phi Angle Controller’ and “PID 
Controller’ are shown in Figures 27 to 29. 
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Figure 27.—Phi_Angle_Controller subsystem. 
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Figure 28.—Trigonometric function Phi Angles. (a) Function 1. (bo) Function 2. (c) Function 3. (d) Function 4. 
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e) Function 5. (f) Function 6. (g) Function 7. (h) Function 8. (i) Function 9. 
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Figure 28.—Concluded. (j) Function 10. 
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Figure 29.—PID Controller subsystem. 


The “Phi Angle” for each function must initially be 
determined at or during actual slow speed (~100 rpm) motor 
operation, in order to define an appropriate value for 
accelerating or decelerating the rotor (either in the clockwise or 
counterclockwise rotation direction). Once values are 
determined, they are then “hardwired” into the code, which will 
automatically apply them during actual slow- or high-speed 
motoring operation. 

Before final implementation of the code, it is probable that 
(depending on the configuration of the motor) a “Phi Angle” 
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scheduling routine is appropriate for improving motoring 
operation as the rotor rotates through progressively higher 
speed values. This scheduling scheme, however, was 
unnecessary in the original code from which the Simulink® 
logic was adopted. In addition, the stator and rotor pole 
combination will dictate the specific range of “Phi Angle” for 
an optimal rate of rotor rotation acceleration or deceleration. 
The following section will discuss the most critical element of 
the control code. 


Wave Shape Generator 


The most critical component of this controller is the 
“Wave Shape Generator” block. It 1s located at the left side of 
the GUI (see Fig. 12). This generator element contains 
10 trigonometric functions that are used to produce the drive 
signals (see Fig. 13) for the motor’s stator coils. The 


“Phi Angle Controller’ and the “PID Controller” are 
simultaneously engaged with the “Wave Shape Generator” 
during normal motor operation and provide the “Wave Shape __ 
Generator” block with the phi angle and “signal amplitude” that 
are used in the argument and amplitude coefficient, respectively, 
in the imbedded trigonometric functions. The “Wave Shape __ 
Generator” subsystem is shown in Figure 30. 


Signal Generator 


theta , aad 
Saturation 
= = 
Nis 5 
Saturation’ 
Frome Ground3 
D es 
Phase3 
ici ) Saturation8 Ground4 
From4 = 
Ground5 


rpm_input 


Saturation4 


Signal Generator 


[C] 
amplitude_input 
Goto4 


Terminator1 
Terminator2 
Terminator3 
Terminator4 


Terminator5 


freq_out 


IS 


function 4 

Saturation6 ~ 

function_5 

Goto1 Saturation7 a. 


function_6 


a(S 
a 
Ww 


nw 
© 
ee 
= 
© 
=: 


IS 


function_7 
2 


” 
- 


uration 


Is 
E 


function_8 
5 


“ 
© 


turation 


@ 
SF 
) 
. 
~ | 
o. 
N 


Multiport 

Switch’ 

=) 
function_select 


Goto2/3 


Figure 30.—Wave_Shape_ Generator subsystem. 
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the left side of the “Wave Shape Generator’ 
(see Fig. 12) there are a number of input ports. The input ports’ 
functions, from top to bottom, are as follows: 


The “theta” port accepts “radians per second values” 
from the “Per Rev Logic” block. 

The “exponent” port accepts constant values for the 
exponent of the Sinu_ Pulse, and Sine functions. 

The “rpm_ input” port accepts the desired input rpm 
values from the “Input (rpm)” block (see top center of 
GUI). 

The “phi” port accepts “Phi Angle Controller’ values 
for phi angle adjustment to the phase angle values. 
The “nl” port accepts values for the maximum number 
of terms in the Fourier series representation for the 
following functions: 

Trapezoidal 01 

Trapezoidal 02 

Square Wave 01 

Square Wave 02 

Sawtooth 

Triangle 02 


Gy, OP, Oe SOT 2 
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e The “function 4” to “function 7” ports accepts values 
to adjust the pulse width of the following functions: 
o  Trapezoidal 01 
o  Trapezoidal_02 
o Square Wave 01 
o Square Wave 02 

e The “function 8” port is unused at this time. 

e The “amplitude input” accepts constant or varying 
amplitude values from the “PID Controller” depending 
on “Sw4” settings. 

e The “function select” accepts values for selecting any 
of the 10 trigonometric function of the 
“Wave Shape Generator.” 


The output ports “Phasel”, “Phase2”, ‘“Phase3”, and 
“freq out” are on the right side of the “Wave Shape _ 
Generator.” Output ports “Phasel” to “Phase3” provide the 
control signals that are ported to the (SPACE hardware which, 
in turn, is used to operate the motor. The motor’s rotor provides 
the pulses that are used in place of “Pulser1” or “Pulser2.” 


S Functions the overall functioning of the controller. Figure 31 shows the 
collection of S Function blocks that are located throughout the 
The green rectangular blocks in the various subsystems are various subsystems. The C codes for the various S Functions 

S Functions. S Functions contain the C codes that are crucial to are listed in the appendix. 
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Figure 31.—Collection of S Functions. 
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Concluding Remarks 


The MathWorks® Simulink® controller logic was adapted 
from an original C controller that successfully drove a room 
temperature, 36-tooth, four-pole, rotor reluctance motor. The C 
logic achieved a required rotor speed of 3000 rpm—the design 
limit of the rotor. The frequency filter in the Simulink® 
controller provided a more accurate input and output signal to 
the rotor by suppressing intrinsic signal anomalies, which were 
determined to be insignificant for high-inertial rotor operation. 
The Simulink® controller is versatile, and with slight 
modifications, can be used to drive other reluctance motor types 
incorporating dissimilar stator and rotor pole combinations. The 
controller was designed with the goal (after appropriate 
modification) of controlling a future superconducting motor 
that will be used as a testbed in the development of other 


NASA/TP—2017-219451 


26 


superconducting aviation propulsion motors envisioned for 
future turboelectric aircraft. The Simulink® findings presented 
in this report only apply in the simulated test environment. 
However, the simulated inputs would be replaced by actual 
proximity probe signals in a real-world dSPACE GmbH 
application. Safe and effective application of this code 
necessitates familiarity with the Simulink® and dSPACE 
technologies. It 1s advised that thorough simulation precede 
actual implementation and application of the code to familiarize 
the operator with its various control options, to prevent any 
unintentional code settings. 


Glenn Research Center 
National Aeronautics and Space Administration 
Cleveland, Ohio, May 17, 2017 


Appendix 


This appendix contains the signal generator functions. The appendix also contains the C programming code for the 
following—signal generator logic, phi angle controller, PID gain selector, radians 1 and 2, and frequency filter 1 and 2. 


Signal Generator Functions 
This section contains signal generator functions. 
1. f(0)=Amp*sin(@) Sine 
2. f(8)=Amp *| 0.5 “(1 + sin(®)) \p,p=9 »Sinu Pulse 


3. f(0) = Amp *| (sin(®))* p |, p =9 > Sine (@)%*p 


4. fs anp*| E14) 2 ys =e 1)" cos(nna) — Jeos( “I 


i > Trapezoidal 01 
2L 


5. poy amo" tO = costa) oos| m=) 


c= > Trapezoidal 02 
2L 


c 2c (-1)' NTC NTUx 
6. x) = Amp*} —-— sin cos} —— || 3% Sq W O01 
f(x) p*—-— (= (== q Wave _ 


, C 
7 = Amp* sin) — meu) in) ——|| > sq Wave 02 
. f(x)= Amp >. sin SE) ED sin I q Wave | 


1 lew 1. 
foo) =Amp*| 2-29" 2 snl 2) epieiiieeth 
1 =n 
1 4 00 ] NIX 
(n- /2 . 
10. f(x) =Amp* Ee =>. ies. —] in| } > Triangle 02 
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Signal Generator Logic 


This section contains the signal generator logic C programming code. 


atic. 1,7 
gouble Amp,2,;82,83,;C2Z;,C3,E1L,b2;,B3, Lb, 7,Cc;p1—3;.141592653539793,ex%1,ex2,1,ex3; 


ai rnc Cio HS Ly 7 2a bw be ck he be eo Sine (8) 
{ 
T= 1/(*rpm/60.0);//-===>. period 
Amp = *amplitude; 
*ul = Amp*sin(*theta (O00 + *pha) * p17 130). 


+ 
*u2 = Amp*sin(*theta + (120 + *phi) * pi/180); 
*u3 = Amp*sin(*theta + (240 + *phi) * pi/180); 
streq out = 1L/T; 
\J//End of 1£(*Function == 1) 


else 
RoW LUNCCLON == 2) 7 ata ceaahisedaa ages Sinu Pulse 
{ 

Z = *exponent; 

T= 1 / (*rpm/60.0);//----> period 


Amp = *amplitude; 
“UL. = Amp pow (0.5 (L.U+simi(* theta - (0 + Ole.) * pa7 leo) 2). 
*u2 = Amp*pow(0.5*(1.0+sin(*theta + (120 + *phi) * pi/180)),z); 
*u3 = Amp*pow(0.5*(1.0+sin(*theta + (240 + *phi) * pi/180)),z); 
*treq out = 1/T; 

\/f(/ene of Lf (*FPunction == 2) 


else 
ie PLU On. =— 2) 77 se ngted nae he Smeanes Sine (0@)%*p 
{ 

Z = *exponent; 

T = 1/ @rpm/ 60.0) 3/7 f=—==> peried 


Amp = *amplitude; 
*ul = Amp*pow(sin(*theta + (000 + *phi) * pi/180),z); 
*y2 Amp*pow(sin(*theta + (120 + *phi) * pi/180),z); 
*u3 = Amp*pow(sin(*theta + (240 + *phi) * pi/180),z); 
*treq out = 1/T; 
\j//Bne Of if (*Function == 3) 
else 
SEF rUuneCeion ==.4) )l lsc eeeewd ea tedeeeaus Trapezoidal 01 
{ 
T = 2*L; 
£ = *rpm/60.0; 
1p = 1717.00) 3 
ec = em 24h; 
Amp = *amplitude; 


n=l; 
exl = OQ; 
ex2 = OQ; 


*Al = pow(n,-2); 
“Bl = pow(=1,;n) *cos (n*pi*cc/ (2*L))=1; 
*Cl = cos(n* (*theta + (000 + *phi) * pi/180)); 
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C2 = cos (n* 
C3 = cos (n* 


El = *Al * *Bi *. *Ci; 
EZ = *AlL * *Bl * CZ; 
Eo = *AlL © 2B *® C33 


exl = 


ex? 
ex3 


n++; 


}//End 


42 = Amo* (0 <5" 
*US = Alp (020 


sex = exl; 
*freq out = pow(2*L,-1); 
}//Bene of if (*Function == 4) 
else 
1 (*funelionm: == S) ili enecccdadnues 


{ 


( 

wii) = Amp* (0.5* (ltec/ (2*1L) ) 
( 
( 


exl + El; 
= ex2 + E2; 
= ex3 + E3; 


7 


of while(n <= *nl1) 


1+cc/ (2*L)) 
1+cec/ (2*L)) 


L = pow(2.0* (*rpm/60),-1); 
CC .= *n*2*L; 


n=O; 
p=0; 
exl = 
ex2 
ex3 = 
Amp = 


0; 
0; 
OF 
* 


amplitude; 


while(n <= *nl1) 


{ 


Nn = 


saa GE 
CZ 
GS 


El 
E2 
K3 


exl 
exZ 
ex3 


pt+; 
}//End 
i = 
*u2 = 
*u3 = 


a Pe 
Laie iy a S| 


}//End o 


(2+ 1) 3 


= pow(n,-2); 
= Cos (n* pi (6C/ (271) ))? 
«(*theta + (000 + *phi) 


( 
= cos(n 
= cos(n* (~theta + (120 + *phi) 
= cos(n 


= *Al * *B1 * *Cl1; 
= *Al * *B1 * C2; 
= *Al * *B1 * C3; 


= exl + El; 

= ex2 + E2; 

= ex3 + E3; 

of while(n <= ¥*nl) 
Amp* (0.5-(4.0/ (pow (pi, 2) * 
Amp* (0.5-(4.0/ (pow(pi,2)* 
Amp* (0'.5=(4,.07 (ow (01,2) 
exl; 


out. = powiZ*b,<—1)+ 
f if (*Function == 5) 
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(*theta + (120 + *phi) 
(*theta + (240 + *phi) 


+2/(pow(pi,2)* 
#2) (DOW (Oiy2) @(1=Ce/ (27). i= e22)-> 
+27 (pow (pL, 2)" (l=ce/ (27) ))* ex3)> 


*(*theta + (240 + *phi) 


* pi/180)) 
~ 94.7 180.) ) 


(l-cc/ (2*L)))* exl); 


ee ew Trapezoidal Oz 


(01/180))); 
(0i/180))); 
* (pi/180))); 


(1.0-cc/L) )) *exl); 
(1 .0=CC/ L):)) 6x2); 
(1.0-cc/L)))*ex3); 
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else 


eee OCC ton == G) 2 aeeseatinnaewadss 


{ 
L = pow(2.0* (*rpm/60),-1); 
Co.= 4m 27h; 
Amp = *amplitude; 


n=1; 

exl = 0; 

ex2 = QO; 

ex3 OF 

while(n <= *nl1) 

{ 
*Al = pow(-1,n)/n; 
“BL = Si7(i* op “cc Lis 
“Cl = Cos (n*(*theta + (000 + *phi) 
C2 = cos (n*(*theta + (120 4+ *phi) 
C3 = cos(n* (*theta + (240 + *phi) 


Hl = *AL = *Bil & *Ci 
E22 = +All * SBl = C2; 
BS = AL * *Bl = Coss 
= exl + El; 
exZ = exZ2 + E2; 
6X3 = 6x3 + Ba? 


OD 

me 

FH 
| 


ay 
}//End of while(n <= *n1) 
*ul = Amp* (cc/L + (2/pi)*exl1); 
s2 = Amp (CC + (2/01) 6x2) > 
*u3 = Amp*(cc/L + (2/pi) *ex3); 


sex = exL; 
*freq out = pow(2*L,-1); 
\/j/enc Gf if (*Function == 6) 
else 
LPC Lon, == 7) Wiens ees esha tees as 


{ 
L = pow(2.0* (*rpm/60.0),-1); 
CO = SZ i; 
Amp = *amplitude; 


n=1; 

exl = 0; 

ex2 = 0; 

ex3 OF 
while(n <= *nl1) 


*Al = pow(-1,n)/n; 


*Bl = Cos (n* 01" cc b)=1; 

“Ci = Sane(n* (“theta + (O00: + “phii) 
C2 = sin(n*~ (*theta + (120 4+ *phi) 
CS = S170 (*theta + (240.4 ~pi) 


BL = *AlL * *Bl * Cl; 
EZ = *Al © *BL * C2; 
Bo = “Al * 4B. Co; 
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..oq Wave O01 


* a) 160) 
* pi 7/180) ) 3 
= a LEO) ) 3 


saanqG Wave U2 


* p1i/180)); 
* p1/180)); 
* p1/180)); 


30 


exl exl + EL; 
ex2 = ex2 + E2; 
ex3 = ex3 + E3; 


De; 
\J//end 6f while(m <= *nl) 
*ul = Amp* (2/pi) *ex1; 
*u2 = Amp* (2/pi) *ex2; 
7S = Amo (27702) *ex3- 


*ex = exl; 
*treq oul = pow(2*h,-1); 
iJ J/End of 2fi* Function == 7 
else 
LE i*funerion ==] OC) le sden dwn sheet aeeae as Sawtooth 


{ 
L = pow(2.0* (*rpm/60),-1); 


) “Ne 


mplitude; 
while(n <= *nl1) 


n= p; 

*Al = pow(n,-1); 

*B1l = sin(n* (*theta + (000 + *phi) * pi/180)); 

B2 san(n*® *thete + (120 + *phi) * pi71s80))-; 
( 


B3 = sin(n* (*theta + (240 + *phi) * pi/180)); 


Bl = *Al * *Bi1; 
F2 = *Al * Ba2; 
BS = *Al * Bo; 


exl exl + EL; 
ex2 ex2 + E2; 
ex3 = ex3 + E3; 


ptt; 
ij /Ene. 6f while(a <= *nl) 
“U1 = Amo (0.5. = pOw(pi,—1) ex); 
sue =. Amp*(U<.0 = DOW (pi, =1)*ex2) > 
“U5 = Amp (Oyo = pow (pi, —1L)*exs) > 


sex = exl; 
*iireg Out = pow(2*4,—-1); 
}//end of Lf (*Function == 9) 
else 
a runner =— 9) 7) 44644. 84044655 68 Sd-d we oS Triangle 01 
t 
Amp = *amplitude; 
i, = pow(2<0* (*rpm/ 60). ,=1) + 


n=OQ; 
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p=0; 


exl = 0; 
ex2 = 0; 
ex3 = 0; 
while(n <= *nl1) 
{ 
= (2 pork) s 
cA = powln;—Z) 7 
*B1l = cos(n* (*theta + (000 + *phi) * pi/180)) 


( 
B2 = cos(m* “theta + (120 + *phi) * pi7180)); 
B3 = cos(n*(*theta + (240 + *phi) * pi/180)) 


El = *Al * *B1; 
B2 = *Al * Ba2; 
BE3 = *Al * = B3; 


exl exl + EI; 
ex2 ex2 + E2; 
ex3 = ex3 + E3; 


ptr; 
\j//End of while(™m <= “nl) 
*ul = Amp* (0.5-(4/pow(pi,2))*ex1) ; 
*u2 = Amp* (0.5-(4/pow(pi,2)) *ex2) ; 
*u3 = Amp* (0.5-(4/pow (pi, 2) ) *ex3); 


*ex = exl; 
*sireg oul. = pow(2*h,;—-L); 
\//End *Function == 9) 
else 
ie (tinction == LO) 7 evnd ac daevndoedw ees Triangle 02 


{ 
L = pow(2.0* (*rpm/60),-1); 


n=0; 

p=0; 

exl. = 0; 

ex2 = 0; 

ex3 = OQ; 

Amp = *amplitude; 


while(n <= *nl1) 
{ 


DPD = (2*pD + -L)g 

*Al = pow(- nie ag cee aa 

*BL = Ssani(n* (*thete + (000 + *phi) * pi7180))-; 
B2 = sin(n*(*theta + (120 + *phi) * pi/180)) 
B3 = sin(n*(*theta + (240 + *phi) * pi/180)) 


EL = SAL * 4B; 
F2 = *Al * Ba2; 
ES. = *AL * Boy 


exl = exl + El; 
ex2 ex2 + E2; 
ex3 = ex3 + E3; 


ray 
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}//End of while(n <= *n1) 

*ul = Amp* (8/pow(pi,2)) *exl1; 
*u2 = Amp* (8/pow(pi,2) ) *ex2; 
*u3 = Amp* (8/pow(pi,2) ) *ex3; 


*ex = exl; 
*freq out = pow(2*L,-1); 
\j/ fend of if (*Fonction =—-10) 


Phi Angle Controller 


This section contains the phi angle controller C programming code. 


[ [KKK KKK KK KKK KKK KK KK KK Acceleration/Braking Block KKKKKKKKKKKKKKKKKKKKKKK 


double speed; 
speed = *rpm2; 


if (speed < *rpmset) 

| *phi = *phi A accel;// Accelerating 

Eeees > *rpmset) 

| *phi = *phi B brake;// Braking 

betenees > (*rpmset- 5)) && (*rpm2 < (*rpmset + 5))) 


{ 
*phi = 0;// Neutral 


} 


[ [KKK KKK KK KKK KKK KK KKK End ‘ona Acceleration/Braking Block KKKKKKKKKKKKKKKKKKKKK 


PID Gain Selector 


This section contains the PID Gain Selector C programming code. 


if (*u0 <= 6000) 
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Radians I and 2 


This section contains the radians | and 2 C programming code. 


‘Lads: = *rads: per _ Sec** cine; 


Frequency Filter 1 and 2 
This section contains the Frequency Filter 1 and 2 C programming code. 
Selector: 
i1f(*p <= 31) 
n=l 


else 
if (*p > .31 && *p <= 62) 


Lf (*ul > 1) 


AVI al 3 
| | ~--------------------- 
u2 > 2 


sa as 
/ | ~--------------------- 
u3 > 3.5 
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